Phase Transitions in liquid Helium 3 
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The phase transitions of liquid 3 He are described by truncations of an exact nonperturbative 
renormalization group equation. The location of the first order transition lines and the jump in the 
order parameter are computed quantitatively. At the triple point we find indications for partially 
universal behaviour. We suggest experiments that could help to determine the effective interactions 
between fermion pairs. 



05.30.Jp, 67.40.Kh, 03.65.Db 

The low temperature phase transitions and the super- 
fluid phases of liquid Helium 3 Jj]] can be described Q by 
a field theory for complex 3x3— matrices A represent- 
ing fermion pairs. An approximation based on quartic 
polynomials in A, including the renormalization group 
running of the coupling constants [||, can account for 
the rough phase structure but not for details of the tran- 
sitions. We will be concerned with the transitions from 
the normal liquid to superfluidity in vanishing external 
magnetic field. In a mean-field treatment those transi- 
tions are of second order, whereas early RG-calculations 
H found hints for them to be fluctuation induced first 
order transitions. The BCS weak-coupling theory esti- 
mates the critical region where fluctuations invalidate the 
mean-field approximation to extend over a temperature 
interval AT w 1CP 8 A". Based on experiments on zero 
sound absorption Q it has been conjectured, however, 
that the critical region might be up to a thousand times 
larger [^). The theoretical answers therefore hinge on 
two problems: first, the "microphysical" effective inter- 
actions at the scale of the fermion pairs (A -1 w 120 A) 
are poorly known and second, the transition from the mi- 
crophysical interactions to the macrophysics of the phase 
transitions (i.e. the thermodynamic potentials) is very 
complex. Only if the second problem can be solved, 
the experimental observations of the macrophysics (phase 
structure, jump in the order parameter, etc.) can be used 
to constrain the microphysical interactions. In this let- 
ter we propose a quantitatively reliable mapping of the 
microphysics onto the properties of the free energy den- 
sity. This is effected by means of extended truncations to 
an exact nonperturbative renormalization group equation 
in contrast to largely uncontrolled approximations 
used earlier [||,fjj. 

The different phases of 3 He are characterized by ex- 
pectation values < A > (order parameters) in differ- 
ent directions: the BW-state, A\ = 1, the ABM-state, 

^2 = yf" (A7 + iX$ + A4 — 1X5), and the planar state, 

A3 = y| i A2 represent the phases of the mean-field dia- 
gram. The Xi are the Gell-Mann matrices. Additionally, 
we examine the direction A4 — »/| (A3 — iXi). Here, 



< A > corresponds to the minimum of the effective po- 
tential U(A) for spatially homogeneous fields. Our main 
aim is therefore the computation of U (A) for a given mi- 
croscopic effective action (at the scale A) for which we 
assume the conventional form j| 

r 4 

T A = / d 3 x { Tr VA^VA + mp + ^ bi h + aK D } 

(1) 

with p = Tr A^ A and Ii the invariants to fourth order al- 
lowed by the symmetry group G = SO(3)xSO(3)xU(l): 

I = P \ h = \Tr A T A\\ h = Tr (A T A)(A T A)* , 
h = Tr {A^A) 2 and h = Tr (A* A){J$A)*. 

We include the dipole interaction 



K D =Tr A* Tr A + Tr {A* A) 



(2) 



as a symmetry breaking perturbation arising from the 
spin-orbit coupling of the atoms in a Cooper pair and 
neglect the strain gradient terms 0. All quantities are 
in units of a characteristic critical temperature T c = 
2.6 mK and m(T) = m(T c ) + + ... reflects the 

temperature dependence. In the paramagnon theory the 
quartic couplings bi are given by pfl 



b = c(2 + .2S), h = -c(l + AS), b 2 = c{2-mS), 

63 = c(2 - .555), 64 = -c(2 + .76) (3) 

with c = 0.001 and 8 parametrizing the pressure depen- 
dence. We study a pressure regime of— 2 < 5 < 1.8. For 
8 > 1.9 the microscopic potential becomes unbounded 
from below. This will be compared to a calculation where 
the couplings are ten times as strong. 

Our approach is based on the effective average action 
Ffe[A(x)] which interpolates between the microphysics 
(k = A) and the macrophysics (k = 0) by means of an 
exact flow equation ||. We truncate the most general 
functional dependence on A(x) by 



r» 



{Z k TrVA^VA + U k {A)+a k K D } (4) 
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where Uk(A) is a function of the eighteen real fields ip y 
in A respecting the symmetry group G. We can compute 
the effective potential U(A) = Uq(A) by a numerical in- 
tegration of the flow equation ||] 

dtU k (A) = ± J 8 t R k Tr ((Z k q 2 + R k ) + Ml(A))- 1 

(5) 

Here, t = In k. The infrared cut-off 

R k (q) = W(l-e^)- 1 (6) 

makes the flow equation ultraviolet and infrared finite. 
We need the mass matrix of second derivatives of U, 
My z = d!p dip • ^ nc now ec l ua tion for the wave func- 
tion renormalization Z k can be found in Q . 

In a first investigation we have used a sixth order poly- 
nomial approximation for U k and searched for minima in 
all the fifteen directions with different residual symme- 
tries given in ^ . In this scheme we found a phase transi- 
tion into the direction A4 for large S. We have convinced 
ourselves, however, that polynomial approximations do 
not allow reliable statements about the phase diagram. 
This is why we deal with an arbitrary dependence oiU k 
on the field in a given direction A = £A n for which we 
write U k (£A n ) = V k (p) with p = 3£ 2 . We assume, that 
the directions n include the absolute minimum of U k . We 
get the flow equation for v£"^ by evaluating (^) at points 
on the line £,A n . To do so, we need to know the second 
derivatives of U k with respect to all fields, however. That 
is why we include in our truncation four coupling func- 
tions b\ (p) multiplying all fourth order invariants other 
than p 2 : 

U(£A n + e) = {VW(p) + J2 tHpM - if V )IU+ £ 

i=l 

(7) 

The numbers are chosen such that the invariants 
Ii — /j p 2 do not contribute on the line £A„. Here e 
denotes fields orthogonal to £A n and we expect that this 
truncation gives a good approximation to the true poten- 
tial in a neighborhood of the line £,A n . 

The flow equations of the coupling functions b[ n \p) 
are derived by forming appropriate linear combinations 
of second and third derivatives of (||) with respect to cer- 
tain fields. At the origin we need to take fourth deriva- 
tives to extract the couplings from our potential trun- 
cation. We have computed these flow equations for the 
four field directions A n . The resulting equations are very 
long, altogether about 4000 lines. They all derive, how- 
ever, from compact expressions, namely field derivatives 
of (J5J) , which can be evaluated by a computer. We con- 
sider this structural simplicity to be one of the major 
strengths of our formalism. 



We have integrated the resulting equations numerically 
by laying the potentials and coupling functions on grids 
of ten points each. In a first run, we detected truncation 
errors in the quantities dv dp \a=q- They originate in an 
insufficient treatment of the symmetry constraints. We 
fight this problem in two steps. First, we infer the first 
potential derivatives at the origin from the term mp and 
integrate a separate flow equation for m. Second, for the 
higher orders we employ for small \t\ a mixed truncation 
of the form 



4 10 

U(£A n + e) = {mp + h h + c j J j + A W (P) 

+ Yp 2 Bl n \pM~ri n) p 2 )}kA n+e (8) 
i=i 

Here Jj are the eleven invariants of sixth order. u/ n ) 
and describe the deviation from a sixth order poly- 
nomial approximation. The symmetry is manifestly re- 
spected up to sixth order. We have estimated the ac- 
curacy of the sixth order polynomial approximation by 
calculating dtU(^A n ) to fourth and sixth order in £ in 
two ways: from the flow of the couplings in the polyno- 
mial part of truncation (|^) and directly by equation (|^). 
The discrepancy stays smaller than 20% as long as the 
flow parameter |i| is not too big, t > —5. It grows to a 
few hundred percent, however, until the phase transition 
occurs at t « —8. This is due to terms of eighth and 
higher order that are truncated in the flow equations for 
the couplings bi and Cj in (||). We conclude, that with 
a polynomial truncation of the free energy to realistic 
order, no conclusive statements can be obtained about 
the phase diagram. As soon as the errors exceed 20%, 
we switch to the truncation (0). The truncation error 
in the quantities dtV^ at the potential minima is esti- 
mated by comparing truncation (Q) with a truncation in 
which we extract the orthogonal masses from the polyno- 
mial truncation in the origin. The errors induced by the 
inaccuracies of dv dp \a=o are of the order of magnitude 
of the truncation errors and both are not bigger than a 
few percent. This truncation should therefore allow us 
to compute the phase structure reliably. 

First, we study the case of large couplings, c = 0.01. 
We neglect the spin-orbit coupling, which is justified in 
this case, the other masses being much bigger. By com- 
puting the temperature dependence of the renormaliza- 
tion of the fourth order couplings, wc find a width of the 
critical region of about 10 -6 -ftT , in good agreement with 
the Ginzburg criterion. The resulting phase diagram is 
shown in figure [D. In a pressure regime of 0.26 < S < 0.46 



2 



a beak of A-phase between the symmetric and the B- 
phase is stabilized by the fluctuations of the order pa- 
rameter. In the middle of this pressure interval the tem- 
perature width of this wedge is about 10 -6 -ftf. 
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FIG. 1. Phase diagram of the effective theory for liquid 3 He 
with strong coupling (c = 0.01). The temperature scale is set 
by AT = 10- 6 K. The same picture is found for c = 0.001 
neglecting the spin-orbit interaction. Then AT = 10 -8 . 

Towards the upper end it rises, however, to lO -5 ^ at 
S = 0.41. All transitions are of the first order. The po- 
tentials at the transitions look qualitatively like figure 



-10 12 
FIG. 2. Pressure dependence of the correlation length at 

the phase transition (in \im, circles) and the discontinuity in 

the NMR-frequency shift (in kHz, squares) 

Figure || shows the pressure dependence of the correla- 
tion length, given by the inverse of the renormalized mass 
at the potential minimum, and the discontinuity of the 
NMR-frequency shift at the phase transition. The latter 
quantity is related to the jump in the order parameter 
squared Ap by 



Av z w 6c • 10 10 Hz ■ Ap 



(9) 



which follows from the mean-field temperature depen- 
dence of Av given in We observe a substantial in- 
crease of the correlation length as we approach the triple 
point S = 0.26. This observation nourishes first hope for 
a fixed point in the vicinity of this point in parameter 



space which would give rise to universal behaviour of the 
system. Accurate predictions could then be made despite 
the uncertainty of our knowledge of the microscopic the- 
ory. Also an examination of the flow of the couplings at 
the origin looks promising for the triple point. We ob- 
serve strong renormalization of the fourth order couplings 
For example, ^j- 4 - w 8, whereas analogous ratios away 
from the triple point, say, at S — 0, grow no larger than 
1.5 at the phase transition . 
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FIG. 3. Flow of the dimensionless renormalized fourth or- 
der couplings near the triple point (c = 0.001, no spin-orbit 
coupling) 

Figure || shows the dimensionless renormalized cou- 
plings bi r . We cannot find a true fixed point behaviour. 
However, ratios of couplings seem to tend to constant 
values, such as ^ = 1, corresponding to a partial fixed 
point. 
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FIG. 4. Phase diagram of the effective theory for liquid 
3 He in the parameter range suggested by the paramagnon 
model. The temperature width of the planar phase is about 
W~ 9 K. 



The qualitative features for smaller couplings c = 0.001 
look similar, with correspondingly smaller AT = 10~ 8 K . 
In this case, however, the dipole term should be included 
(a\ = 2 • 10~ 6 ) which we do next. Figure |] shows the 



3 



phase diagram with first order transition lines. The con- 
siderations in are confirmed by our calculation. The 
width of the stripe of planar phase stabilized by the 
dipole interaction is about 2 ■ 10 _9 -ftT as it was expected 
from mean-field considerations. We have not inspected 
the phase structure much further than 10~ 9 K away from 
the transition line. We expect the same features we found 
above: stabilization of a wedge of A phase with a tem- 
perature width of about 10~ 8 -ftr by the fluctuations. 

Tuning the temperature to the transition from the 
symmetric to the planar phase we find a potential Uk f 
(figure ||) which shows clearly that the transition is of 
first order. (Fluctuations with momenta smaller than kf 
make the potential convex without much influence on the 
equation of state flo|] .) From figure || we can easily infer 
the order parameter jump: Ap « 8- 10~ 5 . This turns out 
to be almost independent of the pressure in the regime 
examined by us. It corresponds to a discontinuity in the 
NMR- frequency shift at the transition of about 50 Hz. 
It was suggested in j| that this should be experimentally 
observable. 
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sition. We show the dependence on £ = y/~£ 
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find Ap w 0.01 for the jump into the B-phase at the triple 
point. This leads to a discontinuity in the frequency shift 
of about 2000iJz as opposed to e.g. 370Hz for c = 0.002. 
We think, that it would be even easier to infer the size of 
the critical region from measuring this order parameter 
jump. At least, this way it should be possible to get up- 
per bounds for the coupling strength that are rather close 
to the range discussed in this article. To facilitate the 
interpretation of such NMR-experiments, we calculated 
the jump of Av for various other values of the coupling 
strength c at 5 — 0.15. We show the result in figure 0. 
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FIG. 6. Dependence of the discontinuity of the frequency 
shift (in kHz) on the coupling strength c 

In conclusion, our quantitative description of the phase 
diagram of 3 He should permit to extract reliable informa- 
tion about the system's microphysical interactions from 
future experiments. 



Let us finally ask how experiments may constrain the 
size of the couplings. We have seen that large couplings 
lead to a widening of the critical region AT. It should 
be possible to exploit this fact by measuring the width 
of the critical region. Varying the temperature at a pres- 
sure right under the triple point, one should be able to 
observe two distinct phase transitions . The width of 
the intermediate phase relates to the coupling strength 
via the Ginzburg criterion. By the described experiment 
one could also measure the pressure corresponding to the 
values S = 0.26 and 6 — 0.46. An approximate form of 
the relation p(S) could be obtained then by linearization. 
Strictly speaking, these arguments hold only if the dipole 
term can be neglected. However, the temperature range 
for the planar phase, 10 -9 -fT, is so narrow, that it would 
probably not even show up in experiments with realistic 
temperature resolution. 

Also the order parameter jump depends strongly on the 
coupling strength. With the larger couplings c = 0.01 we 
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